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Abstract 

The coupled-channels method has been a standard tool in analyz¬ 
ing heavy-ion fusion reactions at energies around the Coulomb barrier. 
We investigate three simplifications usually adopted in the coupled- 
channels calculations. These are i) the exclusion of non-collective 
excitations, ii) the assumption of coordinate independent coupling 
strengths, and iii) the harmonic oscillator approximation for multi¬ 
phonon excitations. In connection to the last point, we propose a novel 
microscopic method based on the beyond-mean-field approach in order 
to take into account the anharmonic effects of collective vibrations. 


1 Introduction 

The field of heavy-ion subbarrier fusion reactions started in the late ’70s, 
when a large enhancement of fusion cross sections was experimentally dis¬ 
covered with respect to the prediction of a simple potential model p]. Even 
though the potential model works well for light systems, such as 14 N+ 12 C, 
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it has turned out that it largely underestimates fusion cross sections for 
heavier systems, such as 16 0+ 154 Sm, at energies below the Coulomb bar¬ 
rier. It has been well recognized by now that this large enhancement of 
subbarrier fusion cross sections is caused by the couplings of the relative 
motion between the colliding nuclei to several nuclear intrinsic degrees of 
freedom, such as low-lying collective excitations in the colliding nuclei as 
well as several nucleon transfer processes mmmm- 

Naturally a standard framework for heavy-ion subbarrier fusion has thus 
been the coupled-channels method [5] by including relevant degrees of free¬ 
dom. This method has not only successfully accounted for the subbarrier 
enhancement of fusion cross sections for many systems but has also provided 
a natural interpretation of the so called fusion barrier distributions HI El- 
In the coupled-channels approach for heavy-ion fusion reactions, the fol¬ 
lowing simplifications are usually employed. Firstly, the model space is 
restricted only to low-lying collective excitations, excluding non-collective 
excitations and giant resonances. Secondly, the coupling strength and the 
excitation energy for each state are taken to be the same as those in an 
isolated nucleus and are assumed to be unaltered during the whole process 
of fusion reaction. Thirdly, a simple harmonic oscillator or a rigid rotor is 
assumed when multiple excitations to higher collective states are involved. 
In this contribution, we shall investigate the validity of each of these as¬ 
sumptions. 

2 Role of non-collective excitations 

Let us start with the first assumption, that is, the role of non-collective 
excitations, which are usually not included in coupled-channels calculations. 

Low-lying collective motions are strongly coupled to the ground state, 
and also have a strong mass number and atomic number dependences. They 
play a major role in heavy-ion subbarrier fusion reactions, and are explicitly 
taken into account in coupled-channels calculations. In addition to the low- 
lying collective excitations, there are many other modes of excitations in 
atomic nuclei. Among them, non-collective excitations couple only weakly 
to the ground state and usually they do not affect in a significant way 
heavy-ion fusion reactions, even though the number of non-collective states 
is large (6j. Couplings to giant resonances are relatively strong due to their 
collective character. However, since their excitation energies are relatively 
large and also are smooth functions of mass number, their effects can be 
effectively incorporated in a choice of internuclear potential through the 
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Figure 1: The quasi-elastic barrier distributions for the 20 Ne+ 90 Zr (the left panel) 
and the 20 Ne+ 92 Zr (the right panel) systems. Here, the quasi-elastic barrier dis¬ 
tribution is defined as D qe \ = —d[da qe \/da^]/dE, where a qe \ and ctr are the quasi¬ 
elastic and the Rutherford cross sections, respectively. These barrier distributions 
are evaluated at the scattering angle of 9\ a b = 150° and are plotted as a function 
of effective energy defined by E e g = 2£^sin(0 c . m ./2)/(l + sin(0 c . m ./2)). The dashed 
lines show the results of the coupled-channels calculations with the collective excita¬ 
tions in the projectile and the target nuclei, while the solid lines take in addition the 
non-collective excitations in the target nuclei into account with a random matrix 
model. The experimental data are taken from Ref. [7j. 


adiabatic potential normalization [lj. 

Although in most of cases, the non-collective excitations do not play a 
significant role unless the system is not very heavy, there are certain cases 
in which they play a role. An example is the 20 Ne + 92 Zr system, for 
which the quasi-elastic scattering at backward angles has been measured 
experimentally [7j. Here, the quasi-elastic scattering refers to the sum of 
elastic, inelastic, and transfer cross sections, and is a counter part of fusion 
reactions [&]. Figure 1 shows a comparison between the quasi-elastic barrier 
distribution, defined as D qe \ = —d[da qe \/dcrR]/dE, where <r qe i and ctr are the 
quasi-elastic and the Rutherford cross sections, respectively, for the 20 Ne + 
90 Zr system and that for the 20 Ne + 92 Zr systems [9] . One striking feature 
is that the experimental quasi-elastic barrier distribution for the 20 Ne + 
92 Zr system is much more smeared than that for the 20 Ne + 90 Zr system 
f?]. The dashed lines in the figure show the results of the coupled-channels 
calculations that include the rotational excitations in 20 Ne as well as the 
collective phonon excitations in 90,92 Zr. This calculation reproduces the 
experimental data for the 20 Ne + 90 Zr system but not for the 20 Ne + 92 Zr 
system. The solid lines, on the other hand, take into account also the non¬ 
collective excitations in 90,92 Zr with a random matrix model [Qj. One can see 




that the smearing of quasi-elastic barrier distribution for the 20 Ne + 92 Zr 
system is now well reproduced by the non-collective excitations of the 92 Zr 
nucleus, whose level density is much larger than that of 90 Zr due to the two 
extra neutrons outside the N = 50 shell closure. 

3 Coordinate dependent coupling strength and 
deep subbarrier fusion hindrance 

The coupled-channels approach expands the total wave function with the 
basis of isolated nuclei. Important inputs for coupled-channels calculations 
are, together with an internuclear potential, the excitation energy and the 
coupling strength for each excitations. Usually the experimental data are 
available for the excitation energy, and the coupling strength can be esti¬ 
mated from a measured electric transition probability [Tj. These values are 
usually employed in coupled-channels calculations assuming that they are 
not altered during the reaction process. This assumption has been examined 
recently by Ichikawa and Matsuyanagi m- They have carried out random- 
phase approximation (RPA) calculations with a two-center shell model po¬ 
tential for the 16 0+ 16 0, 40 Ca-|- 40 Ca, and 16 O+ 208 Pb systems, and have 
demonstrated that the coupling strengths are indeed constant at large dis¬ 
tances but they decrease appreciably in the vicinity of the touching point. 

This implies that the assumption of the constant coupling strength is rea¬ 
sonable for fusion reactions at energies around the Coulomb barrier. How¬ 
ever, at deep subbarrier energies, the inner turning point is close to or even 
inner the touching point mi, and it is important to take into account the 
effect of variation of the coupling strength. Notice that these are the en¬ 
ergies at which the deep subbarrier fusion hindrance phenomenon has been 
observed [3j. In fact, the finding of Ichikawa and Matsuyanagi provides 
a microscopic justification for the damping factor introduced phenomeno¬ 
logically in the adiabatic model for the deep subbarrier fusion hindrance 
phenomenon [12]. 

4 Semi-microscopic modeling of heavy-ion fusion 
reactions with a beyond-mean field method 

In heavy-ion fusion reactions at energies around the Coulomb barrier, multi¬ 
ple excitations to higher members of collective states, such as multi-phonon 
states and high-spin states in the ground state rotational band, often play an 



Figure 2: The results of coupled-channels calculations for the fusion cross sections 
for the 58 Ni+ 58 Ni (the left panel) and the 16 0+ 154 Sm (the right panel) systems. For 
the 58 Ni+ 58 Ni system, the vibrational coupling to multi-quadrupole-phonon states 
are considered in the harmonic oscillator approximation. On the other hand, for 
the 16 0+ 154 Sm system, the rotational coupling within the ground state rotational 
band is taken into account. The experimental data are taken from Refs. mm- 


important role. Figure 2 shows typical examples for this. The left and the 
right panels show fusion cross sections for the 58 Ni+ 58 Ni and the 16 0+ 154 Sm 
systems obtained with the coupled-channels calculations. For the former sys¬ 
tem, the vibrational coupling to quadrupole phonon states are considered, 
while the rotational coupling within the ground state rotational band is 
taken into account for the latter system. One can see that for both the sys¬ 
tems the coupling to the first excited state is insufficient and the couplings 
to the higher members are necessary in order to account for the subbar¬ 
rier enhancement of fusion cross sections. This feature has been demon¬ 
strated beautifully also through the analyses of fusion barrier distributions 

mmm- 

In order to take into account those multiple excitations in coupled- 
channels calculations, one usually uses the rigid rotor model for deformed 
nuclei and the harmonic oscillator model for vibrational nuclei [I]. In reality, 
however, most nuclei have neither a pure harmonic oscillator spectrum nor 
a pure rigid body rotational band, although the rigid rotor approximation is 
reasonable for medium-heavy and heavy deformed nuclei. For example, the 
58 Ni nucleus, which has usually been considered to be a typical vibrational 
nucleus, does not exhibit a level spectrum characteristic to the harmonic 
vibration, e.g., the degeneracy of the two-phonon triplet is considerably 
broken. Moreover, a recent theoretical calculation also indicates that the 
B(E 2) strengths among the collective levels in 58 Ni deviate largely from the 
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Figure 3: The low-lying energy spectra of 58 Ni obtained with the multi-reference 
covariant density functional theory (MR-CDFT) method with the PC-PK1 force. 
The arrows indicate the E2 transition strengths, given in units of e 2 fm 4 . The 
experimental data are taken from Refs. EHG21. 


harmonic oscillator limit EM]. 

There are several theoretical ways to describe an anharmonic vibration. 
Among them, a multi-reference density-functional theory (MR-DFT) has 
been rapidly developed for the past decade tm m [19] . This method is 
based on the so called beyond-mean-field approximation, which incorpo¬ 
rates the angular momentum and particle number projections as well as 
the quantum fluctuation of the mean-field wave function described by the 
generator coordinate method (GCM). 

Figure 3 shows the result of the MR-DFT calculation for the 58 Ni nu¬ 
cleus. The calculation employs the covariant density functional theory 
(CDFT) with PC-PK1 interaction [23] . and we thus call it MR-CDFT. One 
can see that the main feature of the energy spectrum, as well as the E2 
transition strength from 2f to Oj 1 ", are reproduced rather well. It is interest¬ 
ing to notice that the overall pattern of B(E 2) values is quite different from 
what would be expected for a harmonic vibrator, even though the excitation 
energies of the 4+, 2%, and 0 ^“ states are about twice the energy of the 2f 
state. In particular, the E2 transition from the 0^ to the 2^ states is much 
smaller than that from the 4^ and the states to the 2f state. Instead, the 
0 ^ state has a strong transition from the 2 ^ state, which clearly indicates 
that the 0^ state is not a member of the two-phonon triplet. Compared 
to the 0^ state, the E2 transition strength from the O 3 " to the 2f states is 
much larger and is comparable to that from the 4^ and the 2^ states to 
the 2f state. This fact makes the O 3 " state a better candidate for a member 
of the two-phonon triplets, even though the excitation energy is a little bit 












Figure 4: The fusion cross sections (the left panel) and the fusion barrier distri¬ 
butions (the right panel) for the 58 Ni+ 58 Ni system. The dashed line is the result 
of the coupled-channels calculations including the double quadrupole phonon exci¬ 
tations in each 58 Ni nucleus in the harmonic oscillator limit, while the solid line is 
obtained with the multi-reference covariant density functional theory (MR-CDFT) 
method by including the Oj 1 ', 2+, O 3 , 2%, and 4+ states. The dotted line in the left 
panel denotes the result in the absence of the channel couplings. The experimental 
data are taken from Ref. m for the fusion cross sections and from Ref. f!5] for 
the fusion barrier distribution. 


large. 

Notice that, in the harmonic oscillator limit, the B(E2) value from any 
of the two-phonon triplet states to the 2f state is exactly twice the B{E2) 
value from the 2f state to the ground state. The calculated B(E 2) values 
shown in Fig. 3, together with the strong transition from the 2^ to the 
0 ^ states, indicate the presence of large anharmonicity in the quadrupole 
vibrations in 58 Ni. That is, the calculated B(E2) values are significantly 
quenched from the values in the harmonic limit. 

We can now ask ourselves how the deviation of the spectrum from the 
harmonic limit affects the subbarrier fusion reactions of Ni isotopes. Figures 
4 shows the fusion cross section crf us (.E) and the fusion barrier distribution 
D{us{E) = d 2 (Ea ius )/dE 2 [2, 24| for the 58 Ni+ 58 Ni reaction. The dashed 
line shows the result of the coupled-channels calculations including up to 
the double phonon states in the harmonic oscillator limit. All the mutual 
excitations between the projectile and the target nuclei are included. On the 
other hand, the solid line in the figure is obtained with the coupling strengths 
calculated with the MR-CDFT method while adopting a phenomenological 
Woods-Saxon potential. To this end, we include the 0^, 2f, 0^, 2%, and 4^ 
states in the coupled-channels calculations. Again, all the mutual excitation 
channels are taken into account. See Ref. EQj for the details of the calcula¬ 
tions. For a comparison, the figure also shows the result of no-coupling limit 












by the dotted line. One can see that the calculations in the harmonic limit 
overestimate fusion cross sections at the two lowest energies, while the MR- 
CDFT calculations underpredict fusion cross sections around 95 MeV. For 
the energy dependence of fusion cross sections, shown in terms of fusion bar¬ 
rier distribution in the right panel of the figure, the MR-CDFT calculation 
leads to a minor improvement by considerably smearing each peak. 

We have carried out similar calculations also for the 40 Ca + 58 Ni system 
m- We have observed again also for this system that the anharmonicity 
effect in 58 Ni smears the fusion barrier distribution, leading to a better 
agreement with the experimental fusion barrier distribution as compared to 
the results in the harmonic oscillator limit. 

5 Summary 

The coupled-channels approach has been a standard tool in the field of 
heavy-ion subbarrier fusion reactions. We have discussed three assumptions 
used in usual coupled-channels calculations. 

We have first discussed the role of non-collective excitations. Usually, 
coupled-channels calculations take into account a few selected low-lying col¬ 
lective states, neglecting the couplings to non-collective states. We have 
demonstrated that even though the non-collective excitations can indeed be 
neglected in many cases, there are certain systems which show those effects. 
One example is 20 Ne + 92 Zr system, which shows a considerably smeared 
barrier distribution as compared to the 20 Ne + 90 Zr system. By explicitly 
including many non-collective excitations, we have demonstrated that they 
indeed smear the barrier distribution for 20 Ne + 92 Zr system while the effect 
is much smaller for the 20 Ne + 90 Zr system. 

We have next discussed two important input quantities in coupled- 
channels calculations, that is, the coupling strength and the excitation en¬ 
ergy for each state. In most of coupled-channels calculations, these quanti¬ 
ties are assumed to be unchanged during the reaction process and are taken 
to be constants. The recent RPA calculations have indicated that this is the 
case at large distances but they decreases considerably around the touching 
point. We have argued that this leads to important consequences for the 
deep subbarrier hindrance of fusion cross sections. 

We have then proposed a semi-microscopic approach to heavy-ion sub¬ 
barrier fusion reactions. The basic idea of this approach is to combine a 
multi-reference density functional theory (MR-DFT) to a coupled-channels 
calculation. The MR-DFT provides transition strengths among collective 


states without resorting to the harmonic oscillator model or the rigid rotor 
model. The advantages of this approach include i) deviations from the har¬ 
monic limit as well as the rigid rotor limit can be taken into account, ii) 
it can therefore be applied also to transitional nuclei, which show neither 
the vibrational nor the rotational characters, and iii) a natural truncation is 
introduced in the coupling schemes. We have applied this approach to the 
58 Ni+ 58 Ni and 40 Ca+ 58 Ni fusion reactions, and have found that the anhar- 
monicities smear the fusion barrier distributions, somewhat improving the 
agreement with the experimental data. 

One of the important current issues in nuclear reaction theory is to de¬ 
velop a microscopic framework starting from the nucleon degree of freedom. 
It has however been extremely challenging to construct a fully microscopic 
theory which is applicable to heavy-ion subbarrier fusion reactions, thus to 
many-particle quantum tunneling. We believe that the semi-microscopic 
approach presented in this paper provides an important step towards this 
direction. 
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